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ABSTRACT 

We report on the discovery of a new potential galaxy threshing system in the 
COSMOS 2 square degree field using the prime-focus camera, Suprime-Cam, on 
the 8.2 m Subaru Telescope. This system consists of a giant elliptical galaxy 
with My ~ —21.6 and a tidally disrupted satellite galaxy with My ^ —17.7 at 
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a photometric redshift of z « 0.08. This redshift is consistent with the spectro- 
scopic redshift of 0.079 for the giant elliptical galaxy obtained from the Sloan 
Digital Sky Survey (SDSS) archive. The luminosity masses of the two galaxies 
are 3.7 x 1O 12 A^ and 3.1 x 10 9 A / f Q , respectively. The distance between the two 
galaxies is greater than 100 kpc. The two tidal tails emanating from the satellite 
galaxy extend over 150 kpc. This system would be the second well-defined galaxy 
threshing system found so far. 

Subject headings: galaxies: evolution — galaxies: interaction 

1. INTRODUCTION 

Hierarchical clustering scenarios suggest that present day galaxies assembled from much 
smaller building blocks during the course of their evolution (e.g., Peebles 1993). Recent deep 
surveys have found small dwarf galaxies (i.e., building blocks) at z > 2 (e.g., Pascarelle, 
Windhorst, & Keel 1998; Ellis et al. 2001; Taniguchi et al. 2003a; Santos et al. 2004; 
Kneib et al. 2004; for a review see Taniguchi et al. 2003b). At the present epoch, the 
assembly process continues with large galaxies capturing their satellite galaxies. In order to 
understand the mass assembly of galaxies, it is therefore important to investigate such minor 
merger processes in detail. 

It is known that most galaxies in the local universe have satellite galaxies (e.g., Zaritsky 
et al. 1997 and references therein). Satellite galaxies are expected to sink to the center of 
the potential of their host galaxies due to dynamical friction (Ostriker & Tremaine 1975; 
Tremaine 1981); however, satellites with long dynamical friction timescales are still found 
orbiting their hosts. The satellites suffer from tidal disruption as they orbit in the galaxy's 
host dark matter halo. Therefore, the kinematic and morphological properties of the satellites 
can be used to investigate the structure and potential of the host dark matter halo. 

The recent discovery of ultra compact dwarf (UCDs) galaxies in the Fornax cluster 
suggests that such tidal disruption processes can yield small satellite galaxies (Drinkwater 
et al. 2000a, 2000b; Bekki, Couch, & Drinkwater 2001; Bekki et al. 2003). As dwarf 
galaxies orbit around the host galaxy, their outer stellar components are tidally removed. 
This process is also referred to as galaxy threshing. Recently, Forbes et al. (2003) identified 
the first case of a galaxy threshing system in an early release image of the Advanced Camera 
for Surveys (ACS) on the Hubble Space Telescope (Tran et al. 2003; de Grijp et al. 2003). 
The ACS image clearly shows faint, long tidal tails emanating from a satellite galaxy around 
an edge-on spiral galaxy at z = 0.145. 
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Although a pair of tidal tails is often found in major mergers (Arp 1966; Toomre & 
Toomre 1973), it is rare to find such tails in on-going minor mergers. In this respect, 
the discovery of Forbes et al. (2003) provides important information of tidal disruption of 
satellite galaxies around the host galaxy. Following Bekki et al. (2001) we use the term 
"galaxy threshing" for such on-going, tidally-disrupted satellite galaxies. 

Here, we present the serendipitous discovery of a new potential galaxy threshing system 
found in the COSMOS 2-square degree field (Scoville et al. 2006). The most important point 
is that we find clear evidence for the disruption of a small satellite galaxy. Destroyed systems 
are often observed, but faint tidal remnants are generally too faint to be detected in most 
cases. We present the photometric properties of the system and discuss some implications. 
We adopt a flat universe with f2 matter = 0.3, = 0.7, and H = 70 km s -1 Mpc" 1 throughout 
this paper, and use the AB system for optical magnitudes. 

The Cosmic Evolution Survey (COSMOS) is a treasury program on the Hubble Space 
Telescope (HST), awarded a total of 640 HST orbits, to be carried out in two cycles (320 
orbits in cycles 12 and 13 each; Scoville et al. 2006; Koekemoer et al. 2006). This is the 
largest amount of HST time ever, allocated to a single project. COSMOS is a 2 square 
degree imaging survey of an equatorial field in /(F814W) band, using the Advanced Camera 
for Surveys (ACS). 

2. OBSERVATIONS 

The COSMOS HST survey alone cannot be used to address all scientific questions with- 
out additional ground-based and space-based follow-up observations at multiple- wavelength. 
We have been carrying out optical multi-band imaging surveys of the COSMOS 2-square 
degree field centered at a(J2000) = 10 h 00 m 28.6 s and 5(J2000) = +02° 12' 21.0" (Taniguchi 
et al. 2006; see also Capak et al. 2006) using Suprime-Cam, which consists of 5 x 2 CCDs 
with 2k x 4k pixels and has a pixel scale of 0."202 pixel -1 (Miyazaki et al. 2002) on the 
8.2m Subaru Telescope (Kaifu et al. 2000). During our two observing runs in January and 
February 2004, we obtained B, V, r', i', and z' band images of the whole COSMOS field. 

The individual CCD data were reduced and combined using our own data reduction 
software (Yagi et al. 2002) and IRAF. The combined images for individual bands were 
aligned and smoothed with Gaussian kernels to match the atmospheric seeing conditions. 
The FWHM of the PSF of the final images has been matched to 0".92 corresponding to 
the FWHM of the PSF in the z' band image which had the worst seeing of all the data. 
Photometric calibrations are made using spectrophotometric standard stars. Details of the 
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data reduction processes for the COSMOS ACS data are given in Koekemoer et al. (2006). 

3. RESULTS 

In Figure 1, we show our five broad-band images of the galaxy threshing system. The 
threshing system is a pair of interacting galaxies labeled A and B in the B band image (Fig. 
1). Galaxy A appears to be a giant elliptical galaxy while galaxy B appears to be a small 
satellite disk galaxy. It is remarkable that very long tidal tails are emanating from galaxy B. 
These features are quite similar to those found in the first galaxy threshing system reported 
by Forbes et al. (2003). Although we have not yet obtained optical spectroscopy of galaxies 
A and B, their physical association appears to be real, and thus this system is the second 
case of galaxy threshing. In Figure 2, we show an ACS image of this system. In the 

high resolution HST image, there are many other satellite galaxies around galaxy A. This 
group has been identified as ID 10203 in SDSS by Merchan & Zandivarez (2005), They are all 
large, bright galaxies (r' ~ 16 — 17) shown in Figure 3. Note that galaxy B is not contained 
in the SDSS catalog, because the spectroscopic redshift of galaxy B was not obtained. 

3.1. Properties of the Two Galaxies in the Threshing System 

In order to study structural properties of the two galaxies in the threshing system, we 
examine their surface brightness distributions using the IRAF "ellipse" program. In Figure 
4, we show the surface brightness profiles of galaxy A as a function of r 1 / 4 . It follows a de 
Vaucouleur's r 1//4 law, suggesting that it is an elliptical galaxy, consistent with the SDSS 
spectrum. In Figure 5, we show the surface brightness profile of galaxy B as a function of 
radius. It shows an exponential radial profile suggesting that this galaxy is probably a disk 
galaxy. 

Next, we derive the apparent magnitudes of the two galaxies using GALFIT (Peng et 
al. 2002). We apply a de Vaucouleur's r 1 / 4 model for galaxy A and an exponential disk 
model for galaxy B. Since the central region (r < 1") of galaxy A is saturated in our i' band 
image, we have masked out this region. The magnitudes and their errors of both galaxies 
are summarized in Table 1. The scale length of galaxy B and the effective radius of galaxy 
A are shown in Table 2 and Table 3, respectively. 

We also give results for the Sersic fitting in Table 4. The Sersic index of galaxy A ranges 
from ~ 2.5 to ~ 3, being smaller than 4 (which corresponds to a de Vaucouleur profile). Note 
that the central region of galaxy A is saturated in i', and thus the Sersic index may have 
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large uncertainty. On the other hand, galaxy B has Sersic index around 1.4 in each filter. 
These results seem to be basically consistent with our interpretation that galaxy A is an 
elliptical galaxy while galaxy B is a disk galaxy. 

At a redshift of 0.08 (corresponding to a luminosity distance of 363.5 Mpc), we estimate 
an absolute V magnitude (luminosity) of My ~ —21.6 (L v ~ 3.7 x 10 W L V:Q ) for galaxy A, 
and of My 17.7 (L v ~ 1.0 x 10 9 Ly iQ ) for galaxy B. 

3.2. Redshift of the Threshing System 

We found a spectroscopic redshift of z = 0.079 for galaxy A (SDSS J100003.2+020146.4) 
in the SDSS spectroscopic data archive of third data release (Abazajian et al. 2005). Un- 
fortunately, no spectroscopic information on the redshift is available for galaxy B. We de- 
rive a photometric redshift for galaxy B using the \ 2 minimizing method (e.g., Lanzetta, 
Yahil, & Fernandez-Soto 1996). We generate SED (spectral energy distribution) templates 
for < z < 1.0 with a redshift bin of Az = 0.01 using the population synthesis model 
GALAXEV (Bruzual & Chariot 2003). The SEDs of local galaxies are well reproduced by 
models whose star-formation rate declines exponentially (r models); i.e., SFR oc exp(— t/r), 
where t is the age of the galaxy and r is the time scale of star formation. We adopt r = lGyr 
models with a Salpeter initial mass function (a power index of x = 1.35 and a stellar mass 
range of 0.01 < m/M Q < 125) and solar metallicity, Z = 0.02. We then calculate SEDs for 
ages of t — 1, 2, 3, 4, and 8 Gyr, corresponding to the SED of a starburst, Irr, Scd, Sbc, 
and elliptical galaxy, respectively. We adopt a visual extinction of Ay = 0. Note that our 
threshing system consists of basically nearby galaxies. Therefore, our results should not be 
affected by metallicity and Ay. Applying this method, we obtain a photometric redshift of 
the galaxy B as z^ = 0.08, being consistent with the spectroscopic redshift of the galaxy A. 
For the rest of the paper, we adopt z = 0.08 for the galaxy threshing system (1" ~ 1.76 kpc 
at z = 0.08). 

3.3. Properties of the Tidal Tails 

We investigate the morphological properties of the two tidal tails emanating from galaxy 
B. As shown in Figure 1, the two long tidal tails are obvious in images of all bands, suggesting 
that they consist of stellar material. The lengths of the two tidal tails are 58 arcsec for the 
tail extending towards the galaxy A and 30 arcsec for the counter tail (east of galaxy B). 
These values correspond to lengths of 102 kpc and 52 kpc, respectively in the plane of the 
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sky. 

The tail extending toward galaxy A appears to be narrower than the counter tail whose 
width increases with increasing distance from galaxy B. This property is also seen in the 
first galaxy threshing system found by Forbes et al. (2003). In order to study this tendency 
more quantitatively, we examine the surface brightness distributions of the tails. In Figure 6, 
we show the i'-band surface brightness profile perpendicular to the tails for the four spatial 
positions, a, b, c, and d that are indicated in the right panel of Figure 6. As shown in the left 
panel of Figure 6, the tail extending to galaxy A is significantly narrower than the counter 
one. The tail widths measured at 3a and 5a noise level for each position are given in Table 
5. The width at position d exceeds 10 kpc even then it is measured at the 5a noise level. 
However, the width of the tail extending to galaxy A is only ~ 6 - 7 kpc. 

There are some knotty structures in the tidal tails. It is likely that some of these may 
be tidal dwarf candidates (e.g., Zwicky 1956; Schweizer 1978; Barnes & Hernquist 1992; 
Yoshida, Taniguchi, & Murayama 1994; Due & Mirabel 1998 and references therein). Future 
spectroscopic identification will be necessary to confirm them. 

We estimate the mean surface brightness of the tidal tails in a sky area above 1.5a of 
the surface brightness in each band. The errors are estimated using following relation; 

(mean counts/pixel) + ct2 ack d 
a (mean counts/pixel) = — - . 

iV is the total number of pixels of tidal tail, and ^background is la sky background noise. The 
results are given in Table 6. The mean surface brightness is fainter by 1 magnitude than 
that of the tail found by Forbes et al. (2003). 

We also estimate the total magnitude of the tidal tails in each band. Errors are estimated 
using following relation; 

a 2 (total counts) = Total counts + N x o"b ackground . 

The results are given in Table 6. It is found that the magnitude of the tails is comparable 
to that of galaxy B in each band. 

In Figure 7, we show the B — r' color distribution in the galaxy threshing system. The 
mean B — r' color of galaxy A is 1.37 ± 0.06, while, that of galaxy B is 1.02 ± 0.07. There 
is a small color difference between the main body of galaxy B and the two tidal tails. This 
suggests that stars located in the tidal tails were dispersed from the main body of galaxy 
B, providing strong evidence that galaxy threshing has been occurring in this system. The 
total magnitude of the tidal tails is comparable to that of the main body of galaxy B. We 
estimated that the luminosity of galaxy B is L v ~ (1.03 ± 9.48 x 10~ 5 ) x 10 9 Ly iO , and that of 
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the tail is L v ~ (0.94 ±2. 16 x 1CT 3 ) x lO 9 Ly j0 . These luminosities are nearly the same. This 
suggests that galaxy B lost about half of its former mass during the (on-going) threshing 
event. 



We assume that galaxy A has a mass-to- light ratio of 10 typical for an elliptical galaxies 
and 3 for galaxy B. Given the V^-band luminosity in section 3.1, we obtain Ma ~ 3.7x lO 12 ^© 
and M B ~ 3.1 x 1O 9 A^ for A and B, respectively. 

Next, we estimate the probable pericenter distance (r p ) for the galaxy threshing event. 
As shown in Figure 5, galaxy B shows an exponential surface brightness profile for r < 2.9 
arcsec. If we assume that the mass outside of this radius has been expelled during the 
previous encounter between the two galaxies, and we use the equation of tidal radius as a 
rough guide of the pericenter distance, we can then estimate the pericenter as follows, 



Using the derived masses of the two galaxies, we obtain a pericenter distance of r p ~ 68 
kpc. Bekki et al. (2003) made general numerical simulations of a galaxy threshing systems 
with r p = 65 kpc. Mayer et al. (2001) used a similar value, r p = 75 kpc. Our estimate is 
consistent with the results of these numerical simulations. 

The projected separation between two galaxies is 102 kpc. If we adopt an average 
relative velocity between the two galaxies of 200 km s _1 , we find that the previous encounter 
occurred r ~ 5 x 10 8 (v / 200 km s _1 ) yrs ago. 

Galaxy B found in this study will merge into galaxy A after a few billion years after 
(see Bekki et al. 2001). Even in a facial phase of such minor mergers, faint remnants could 
be found around the host galaxy. 

Recently, morphological evidence for such a final phase of minor mergers has been 
found around field red elliptical galaxies (van Dokkum 2005). Deep optical imaging survey 
for ordinary- looking galaxies should provide more information on the global, and dynamical 
evolution of galaxies. 



3.4. Dynamical Properties of the Galaxy Threshing System 
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4. Concluding Remarks 

In this paper, we have presented the discovery of a new potential galaxy threshing 
system. Although merging systems and interacting systems are often observed, it is rare to 
find faint tidal tails in on-going minor mergers because their tidal tails are generally too faint 
to be detected. Our discovery is important because we find clear evidence for the disruption 
of a small satellite galaxy. This is only the second well-defined galaxy threshing system 
found so far. 
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Table 1. Basic properties of the threshing system. Magnitudes and their errors are 

derived using GALFIT. 



Name 


RA(J2000) 


DEC(J2000) 


B 


V 


r' 






J095959.6+020206(B) 


9 59 59.61 


+2 02 06.6 


20.87 


20.15 


19.85 


19.56 


19.28 








±0.01 


±0.01 


±0.01 


±0.01 


±0.01 


J100003.2+020146(A) 


10 00 03.23 


+2 01 46.5 


17.29 


16.33 


15.92 


15.48 


15.21 








±0.00 


±0.00 


±0.00 


±0.00 


±0.00 



Table 2. Scale length of the threshed galaxy (galaxy B). These values are derived using 

GALFIT. 



Name 




B 


V 


r' 


i' 


z' 


J095959.6±020206 


-R s (arcsec) 
-Rs(kpc) 


0.50 
0.88 


0.58 
1.02 


0.56 
0.99 


0.53 
0.93 


0.57 
1.00 



Table 3. Effective radius of the threshing galaxy (galaxy A). These values were derived 

using GALFIT. 



Name 




B 


V 


r' 


%' 


z' 


J100003.2±020146 


i?e(arcsec) 
-Re(kpc) 


3.29 
5.97 


3.00 
5.29 


3.20 
5.64 


2.34 
4.12 


4.12 
5.16 



Table 4. Sersic parameter of both galaxies. These values are derived using GALFIT. Note 
that the central region (r < 1 arcsec) of galaxy A is saturated in %' . 



Name B V r' i' z' 



galaxy A 2.91 2.56 2.69 1.80 2.57 
galaxy B 1.41 1.33 1.46 1.47 1.50 
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Table 5. Width of the tidal tails. The positions of the cuts across the tidal tails are 

indicated in Fig.6. 







a 


b 


c 


d 


3a 


(kpc) 
(arcsec) 


7.6±8:i 

4.31°;° 


7 q+0-3 

4 -°-0.2 


11.8l|i 

6.718:? 


i2.3i8i 
7.ol8:l 


5(7 


(kpc) 
(arcsec) 


3.618:8 


c c;+0.2 
°- d -0.1 

3 1+0-1 


8.5i8i 

4.8l8i 


io.3l8:2 
5-918:1 



Table 6. Properties of the tidal tails. 



B V 



Mean surface brightness (mag/arcsec 2 ) 27.90 ± o.os 27.30 ± 0.04 27.11 ±0.07 26.84 ± 0.04 26.29 ± 0.0s 
Total magnitude (mag) 20.91 ± 0.07 20.18 ±0.05 19.82 ± o.os 19.66 ± o.o 4 19.23 ± o.oe 

Background (mag/arcsec 2 ) 29.89 29.39 29.39 28.98 28.05 
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Fig. 1 Suprime-Cam images of the threshing system. Image size of these figures are 2.3' x 1.5' 
(corresponds to 243 kpc x 158 kpc at z ~ 0.08). 





Fig. 2 HST ACS / — 814 band image of the galaxy threshing system. Image size of this 
figure is 2.3' x 1.5'. 
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Fig. 3 Circled galaxies are the member of the galaxy group. They were identified by 
Merchan & Zandivarez (2005). 
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Fig. 4 Surface brightness profiles of the threshing galaxy. Note that the central region 
(< 1") of the threshing galaxy is saturated in %' and has been masked out. We also show the 
la of sky background and seeing radius (= 0.46"). The values of background are shown in 
the Table 6. 
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Fig. 5 Surface brightness profiles of the threshed galaxy. We also show the la of sky 
background and seeing radius (= 0.46"). The values of background are shown in the Table 
6. 
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Fig. 6 Surface brightness distribution perpendicular to the extent of the tidal tails as a 
function of position along the tails. The investigated spatial positions are marked by bars in 
the right panel. The surface brightness profile at each position is shown in the left panel. 
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Fig. 7 B — r' color map. Image size of this figure is 2.3' x 1.5'. Outer parts of the threshing 
system (in particular, galaxy A) show redder colors. However, this feature is due to poor 
signal-to-noise ratio because of much lower surface brightness. 
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